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Abstract 

A detailed analytic and numerical study of baryogenesis through leptogenesis is per- 
formed in the framework of the standard model of electroweak interactions extended 
by the addition of three right-handed neutrinos, leading to the seesaw mechanism. 
We analyze the connection between GUT-motivated relations for the quark and 
lepton mass matrices and the possibility of obtaining a viable leptogenesis scenario. 
In particular, we analyze whether the constraints imposed by SO(10) GUTs can be 
compatible with all the available solar, atmospheric and reactor neutrino data and, 
simultaneously, be capable of producing the required baryon asymmetry via the lep- 
togenesis mechanism. It is found that the Just-So 2 and SMA solar solutions lead to 
a viable leptogenesis even for the simplest SO(10) GUT, while the LMA, LOW and 
VO solar solutions would require a different hierarchy for the Dirac neutrino masses 
in order to generate the observed baryon asymmetry. Some implications on CP 
violation at low energies and on neutrinoless double beta decay are also considered. 



1 Introduction 



Baryogenesis via leptogenesis is one of the most appealing mechanisms for 
generating the observed baryon asymmetry of the Universe (BAU) . The inter- 
est in leptogenesis has been reinforced after the strong evidence for neutrino 
oscillations reported by the Super-Kamiokande experiment [1] and recently 
confirmed by the results of the Sudbury Neutrino Observatory (SNO) [2], 
both pointing towards nonzero neutrino masses. From a theoretical point of 
view, the smallness of neutrino masses can be naturally explained through the 
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seesaw mechanism [3] which automatically follows, once right-handed neutri- 
nos are added to the standard model (SM). Since the right-handed neutrino 
fields are singlets under the SU(3) C x SU(2) x U(l) gauge symmetry of the 
SM, right-handed neutrino Majorana masses are not protected by the gauge 
symmetry so they can be much larger than the electroweak scale. These right- 
handed neutrino states appear naturally in Grand Unified Theories (GUT) 
such as SO (10) [4], where they are included in the same irreducible represen- 
tation together with quarks and leptons. 

At low energies, all the information about neutrino masses and mixing is 
contained in the effective neutrino mass matrix which can be expressed by the 
seesaw formula: 

M v ~ -M D M- R X M T D , (1) 

where Mo and Mr are the Dirac and right-handed neutrino mass matrices, 
respectively. Unfortunately, M D and M R are not fully determined by the low- 
energy experimental data. In this respect, GUTs can play an essential role 
since in their framework Mo is closely related to the quark or charged lepton 
mass matrices. Another alternative to this approach is to impose some extra 
symmetries like horizontal or discrete symmetries [5]. 

In this paper we investigate the constraints on Mo resulting from the require- 
ment of a viable leptogenesis [6-8]. More specifically, we investigate whether 
leptogenesis can discriminate among the various solar neutrino solutions in the 
framework of some GUT- inspired patterns for M D . We shall assume a special 
form for Mo in the weak basis where Mr and the charged leptonic mass ma- 
trix Mi are diagonal. There is of course no loss of generality in choosing these 
two matrices real, positive and diagonal. Since our discussion does not rely on 
specific textures for these matrices it will cover a wide range of mass matrices. 
In this sense our analysis is quite general. Yet, as we shall see, the hierarchy 
of the masses resulting from Mo plays a crucial role in the viability of the 
leptogenesis scenario. In fact, mass hierarchies such as the ones imposed by 
minimal SO(10) strongly constrain the allowed solar solutions. Nevertheless, 
it turns out that allowing for a more general choice of the mass spectrum of 
Mo, it is possible to reconcile the different solar solutions with the required 
cosmo logical BAU. 

The connection between neutrino masses and mixing and the observed baryon 
asymmetry of the Universe becomes clear if one takes into account that the 
leptogenesis mechanism starts with the CP asymmetry generated through 
the out-of-equilibrium L-violating decays of the heavy Majorana neutrinos 
[9,10,6], leading to a lepton asymmetry which is subsequently transformed into 
a baryon asymmetry by the (B + L)-violating sphaleron processes [11]. As it 
will be discussed later, this baryon asymmetry depends mainly on the heavy 
Majorana neutrino mass spectrum and the Dirac neutrino Yukawa coupling 
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matrix, in the basis where the mass matrix Mr is diagonal. These matrices 
will be constrained by GUT-like relations and by the available information on 
neutrino mixing and mass squared differences measured in solar, reactor and 
atmospheric neutrino experiments. On the other hand, requiring the produced 
baryon-to-entropy ratio Yb = ub/s to be in the presently allowed range [12], 

1.7 x KT 11 <Y B <8.1 x KT 11 , (2) 

can in principle constrain some mixing and mass parameters, or even discrim- 
inate the solar neutrino solutions which are compatible with the leptogenesis 
scenario. We shall see that the only parameters relevant to this analysis are the 
absolute value of the (l,3)-element of the leptonic mixing matrix, the Dirac- 
type phase and the two Majorana phases appearing in this matrix, as well as 
the mass of the lightest neutrino. 

At this stage one may wonder whether there is any link between leptogenesis 
and leptonic CP violation at low energies since they both arise from the 
phases appearing in the leptonic mixing matrix. Several authors have already 
addressed this question under different assumptions [13,14]. Experimentally, 
leptonic CP- violating effects can be probed in neutrino oscillation experiments 
and they are only sensitive to the Dirac-type phase appearing in the leptonic 
mixing matrix. On the other hand, the size of double beta decays is affected 
by the presence of Majorana- type CP- violating phases. 

The paper is organized as follows. In Section 2 we present the general frame- 
work and the strategy to follow. Special emphasis is given to the construction 
of the Dirac, Majorana and effective neutrino mass matrices as well as to the 
identification of the CP- violating phases relevant to leptogenesis. In Section 3 
we present the simplest leptogenesis scenario, namely, the out-of-equilibrium 
decay of a single heavy Majorana neutrino. We briefly recall how to estimate 
the baryon asymmetry taking into account the washout effects due to inverse 
decays and lepton number violating scattering processes. This section is di- 
vided into two parts. First we develop a simple analytic approach which will 
allow us to obtain upper bounds on the baryon asymmetry for the different 
solar neutrino solutions and assuming different patterns for the light neutrino 
spectrum. We then perform a detailed numerical analysis in order to identify 
the regions of the parameter space where the leptogenesis mechanism is effi- 
cient enough to produce the required cosmological baryon asymmetry, taking 
into account the presently available low energy neutrino data. Section 4 is 
devoted to discuss some implications on leptonic CP violation at low energies 
and neutrinoless double beta decay. Finally, our conclusions are presented in 
Section 5. 
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2 General framework 



We shall work in the framework of a minimal extension of the SM with only one 
right-handed Majorana neutrino per generation. After spontaneous symmetry 
breakdown, the leptonic mass terms are of the form 

-C mass = 11 M e £ ° + v\ M D u R + X - u R T CM R u R + h.c. , (3) 

where £ ®,r an d v \r are the charged lepton and neutrino weak eigenstates, 
respectively. The charged lepton mass matrix Mg and the Dirac neutrino mass 
matrix are in general 3x3 complex matrices whilst the heavy Majo- 
rana neutrino mass matrix Mr is constrained to be symmetric. The right- 
handed Majorana mass term is SU(2) xU(l) invariant, as a result it is the only 
fermionic mass term present in the Lagrangian before the symmetry breaking. 
The Lagrangian in Eq. (3) can be rewritten as 



-Cr. 



niCM*n L + I°M e e R , + h.c, 



(4) 



with n L = (y° L , and 



M 



Ml 




(5) 



The zero entry in A4* is due to the absence, at tree level, of a Majorana mass 
term for the left-handed neutrinos v L . The diagonalization of M* is performed 
through the transformation 



V T M* V = V = diag( m 1 ,m 2 , m 3 , M u M 2 , M 3 ) , 



(6) 



where and M; denote the physical masses of the light and heavy Majorana 
neutrinos, respectively. It is useful to write V and D in a block form as 



V = 



V 



T 



V = 



d u 



\ 



d 



(7) 



Rj 



with d v = diag(mi, m 2 , m 3 ), dn = diag(Mi, M 2 , M 3 ). From Eq. (6) one derives 
the approximate relations 



-K^MdM^MIK* = K^M V K* = d v , 



(8) 



S-t = -K^M D M R X , 



(9) 
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together with the exact relation 



R = M D T*d R 1 . (10) 

The relation in Eq. (8) is the usual seesaw formula, which leads in a natural 
way to small left-handed Majorana neutrino masses of order M D M R l . The 
matrix K is almost unitary, which is in agreement with the matrices S and R 
being of order Md M R l as can be inferred from Eqs.(9) and (10). 

The neutrino weak eigenstates z/° L are related to the mass eigenstates by 

/ 



Va, = V ia v aL = (K,R) 



N lL 



i = l,2,3, a = l,2,---,6, (11) 



and thus the leptonic charged current interactions are given by 

(h !,Ku L + l L 7m R N L ) W» + h.c. , (12) 

so that K and R give the charged current couplings of the charged leptons to 
the light neutrinos I/; and the heavy ones N i: respectively. Without imposing 
any restriction on the mass matrices, Eq. (3) can be rewritten in the weak 
basis (WB) where M e and M R are real and diagonal in the form 

-C mass = II d e £ ° + v\ M D v\ + l - u° R T Cd R u° R + h.c. , (13) 

where de = diag(m e , m M , m T ). We have kept for simplicity the same notation 
of Eq. (3) for the rotated fields and mass matrices in the Lagrangian (13). 
Obviously, in this WB all CP- violating phases appear in Mp. Since for n 
generations there are n(n — 1) independent phases [14,15], the matrix Mjy has 
six CP-violating phases. In addition, the Lagrangian C mass has fifteen real 
physical parameters: nine contained in Md plus the six masses in de and d R . 

It is easy to see how these phases appear in Mb- Using the polar decom- 
position, Md can be written as the product of a unitary matrix U times a 
Hermitian matrix H. Factoring out all possible phases we can write 

M D = UH = P ( U'P 1 H'P 2 , (14) 

with P^ = diag(e 1 ^ 1 , e^ 2 , e^ 3 ), Pj = diag(l, e ltp i, e l<p z) and the matrices U' and 
H' containing only one phase each. Since P^ can always be rotated away by 
a simultaneous phase transformation of the left-handed charged lepton fields 
and left-handed neutrino fields, only six phases are physically meaningful. 

In the physical basis all CP- violating phases are shifted to the leptonic mixing 
matrix and thus will appear in the matrices K and R. Since to a very good 
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approximation K is unitary and three of its phases can be rotated away, we end 
up with three physical phases in K, two of which are of Majorana character. 
These phases are in general complicated functions of the six phases of Md and 
can be obtained from Eq. (8) by replacing M R by d R . The matrix R, on the 
other hand, verifies Eq. (10). In the WB where M R is real and diagonal, the 
matrix T in (7) is close to the identity. As a result, from Eqs. (10) and (14) 
we obtain 

R = U'P l H'd R l P 2 . (15) 

Once again two of the phases appearing in R are of Majorana character - the 
ones contained in P 2 . As we shall see in the next section, leptogenesis is only 
sensitive to the phases appearing in M D M D , or equivalently R^R. Therefore, 
from Eq. (15) we can conclude that the only CP-violating phases relevant for 
leptogenesis are the phase contained in H' and the two Majorana phases of 

In the exact decoupling limit, R can be neglected and only K is relevant. In 
this limit, the leptonic mixing matrix K is usually referred as the Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) matrix [16,17], which from now on we shall 
denote as U v . However, since we are interested in studying the connection 
between CP violation at low energies and leptogenesis we should keep both 
K and R or, alternatively, the Dirac neutrino mass matrix M D in the WB 
where Mi and M R are diagonal. 

Experimentally, only the charged lepton mass spectrum is very well known. In 
contrast, the experimental data on light neutrino masses and leptonic mixing 
define different allowed regions in parameter space, still leaving a lot of freedom 
in the choice of patterns for the matrices U v and M v . The heavy neutrino 
masses and their mixing with the charged leptons are even less constrained. 

From the theoretical point of view, grand unified theories such as SO (10) 
[4] are a suitable framework not only to analyze fermion masses but also to 
implement the seesaw mechanism. One of the attractive features of the SO (10) 
model is that its gauge group is left-right symmetric and, consequently, there 
exists a complete quark-lepton symmetry in the spectrum. In particular, the 
fact that all left-handed (right-handed) fermions of each family fit into the 
single irreducible spinor representation 16 (16) of SO(10) and that the right- 
handed neutrino is precisely contained in this representation is remarkable. 
Several constraints on fermion masses are usually implied in these models 
[4]. For instance, if there is only one 10 Higgs multiplet responsible for the 
masses, then we have the relation M u = M d \ I, M D (the indices u and 
d stand for the up and down quarks, respectively). Similarly, the existence 
of two 10 Higgs multiplets implies Mi = M d and Md = M u . On the other 
hand, if the fermion masses are generated by a VEV of the 126 of SO(10), 
then the SU(4) symmetry yields the relations 3 M d = —Mi and 3 M u = —M D . 
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Of course, equalities such as the ones arising purely from the 10-dimensional 
representation cannot be exact in realistic models, since they imply unphysical 
relations among the quark and charged lepton masses. Additional assumptions 
are therefore necessary in order to predict the correct fermion spectrum [18- 
20]. 

Without loss of generality, one can choose a WB where and Mr are diag- 
onal and real while Md is an arbitrary complex matrix which we write as: 

M D = Vld D U R , (16) 

where dr> = dia.g(m D i,m D 2,m D3 ) and Vl, Ur are unitary matrices. In a min- 
imal SO (10) scenario it is expected a small misalignment between the Her- 
mitian matrices MgM[ and MpM^, similar to that in the quark sector. The 
matrix Vl is analogous to the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
of the quark sector and under the above assumption of small misalignment, 
it should be close to the identity matrix. The Dirac neutrino mass spectrum 
will be constrained by SO(10) relations such as the ones discussed above. For 
defmiteness, we will assume that dr> is given by the up-quark spectrum, i.e 

d D = diag(m n , m c , m t ) . (17) 

A brief discussion of the most general case with arbitrary masses is presented 
at the end of Section 3.2. 

In our analysis we choose d v and U v in agreement with the present exper- 
imental data. The corresponding effective neutrino matrix M u can be then 
computed from Eq. (8), i.e. 

M v = XJ v d v U T v . (18) 

We shall consider different patterns for the PMNS matrix U u corresponding to 
the different solar solutions (LMA, SMA, LOW, VO, Just-So 2 ). For each solu- 
tion we impose the corresponding mass squared differences Am^ = Am 2 2 = 
\m\ — m 2 | (solar neutrinos) and Am 2 = Am| 3 = |m| — m|| (atmospheric 
neutrinos) to be in the presently allowed experimental range. We also let the 
mass of the lightest neutrino to vary in the appropriate range: its lower bound 
is typically fixed by requiring the heaviest Majorana mass to be below the 
Planck scale, while the upper bound is chosen so that neutrino less double 
beta decays experiments are not violated. 

An important feature of this approach is that once M v is chosen, and for a 
givenP] Vl , the masses of the heavy Majorana neutrinos are fixed together 



2 Throughout our paper we either take Vl equal to the identity or close to it. 
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with the matrix Ur, which appears in Eq. (16). Indeed, the relations 



M v = -M D d~ R l M D T = -Vl d D Md D V[ , 
M = -dn 1 V L M u V L T d D l = U R d^ U R T , (19) 

imply that the matrix M is fully determined. Obviously, we have M M* = 
Ur d^ 2 U R , corresponding to an eigenvalue equation. This equation determines 
the heavy Majorana spectrum dR and also the unitary matrix Ur up to a phase 
ambiguity, which can be then fixed by going back to Eqs. (19). 

Notice that Eqs. (19) clearly illustrate a remarkable feature of these kind of 
models [21], which is the fact the large mixing required for the atmospheric 
solution as well as some of the solar solutions can be obtained with small 
mixing in Ur, so that nearly maximal neutrino mixing is produced through 
the seesaw mechanism. In other words, starting from nearly diagonal mass 
matrices and as a result of the interplay between the small off-diagonal entries 
of Ur and a strong hierarchy in do and dR, it is possible to obtain a large 
neutrino mixing via the seesaw mechanism. 



3 Heavy Majorana neutrino decays, leptogenesis and neutrino os- 
cillations 



The crucial ingredient in leptogenesis scenarios is the CP asymmetry gener- 
ated through the interference between tree-level and one-loop heavy Majorana 
neutrino decay diagrams. In the simplest extension of the SM, such diagrams 
correspond to the decay of the Majorana neutrino into a lepton and a Higgs 
boson. Considering the decay of one heavy Majorana neutrino iVj, the CP 
asymmetry in the SM is then given by 

_ T(N i -+lH)-T(N i -+lH*) 
Nl T(Ni^lH)+T (Ni -> IH*) ' 1 j 



If the heavy Majorana neutrino masses are such that Mi < M2 < M3 only 
the decay of the lightest Majorana neutrino N\ is relevant for the lepton 
asymmetry^ and one obtains [23,6] 



• TTV 



D)U i=2,3 



£ Im \{M^M D ) 2 H 



'ill 2 



H4 + ° w, 



(21) 



3 This is a reasonable assumption if the interactions of the lightest Majorana neu- 
trino N% are in thermal equilibrium at the time of the decays, so that the 
asymmetries produced by the heaviest neutrino decays are erased before the light- 
est one decays, or if are too heavy to be produced after inflation [22]. 
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where is the Dirac neutrino mass matrix in the basis where Mr is diagonal 
and v = (H )/y/2 ~ 174 GeV. The funct ions f(x) and g(x) denote the one- 
loop vertex and self-energy corrections, respectively, and are given by 



/(*) 



1 + (1 + x) In 



+ x 



g{x) 



y/x 



X 



(22) 



In the limit Mi <C M%, M3 the CP asymmetry (21) is approximately given by 



3 / M, Mr 
16 nv 2 \ M 2 M 3 / 



(23) 



where 



/ii = 



Im 



(M^Md)^ 



(24) 



(Mj, M D )n 

The lepton asymmetry Y L is related to the CP asymmetry through the relation 



v n L -n L e Nl 
Y L = = d — , 

s 9* 



(25) 



where g* is the effective number of relativistic degrees of freedom contributing 
to the entropy and d is the so-called dilution factor which accounts for the 
washout processes (inverse decay and lepton number violating scattering). In 
the SM case, g* = 106.75. 

The produced lepton asymmetry Yl is converted into a net baryon asymmetry 
Y B through the (B + L)-violating sphaleron processes. One finds the relation 
[24,6] 



Y B = £Y l 



^-1 



Yl , i 



8N f + 4N H 
22 N f + 13 N H 



(26) 



where Nf and Nh are the number of fermion families and complex Higgs 
doublets, respectively. Taking into account that Nf = 3 and N H = 1 for the 
SM, we get £ ~ 1/3 and 



(27) 



The determination of the dilution factor involves the integration of the full 
set of Boltzmann equations. A simple approximated solution which has been 
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frequently used is given by [25] 



v / (U^exp(-|^OTK) , k>10 6 
0.24(k hiK)- 3 / 5 , 10<k<10 6 
l/(2«) , 1<k<10 



(28) 



0<*<1 



where the parameter k, which measures the efficiency in producing the asym- 
metry, is defined as the ratio of the thermal average of the Aq decay rate and 
the Hubble parameter at the temperature T = Mi , 



M f 



K 



(M^M D ) U 



1.7 x 8nv< 



'9* 



Mi 



(29) 



Mp ~ 1.22 x 10 19 GeV is the Planck mass0. A slightly modified approximation 
is used instead by the authors of Ref . [26] : 



0.30(k hiK)- 3 / 5 



10 < k < 10 6 



1/(2 V / ^T9) , 0<^<10 



(30) 



However, it has been recently pointed out [8] that in some cases the above 
approximations could seriously underestimate the suppression in the baryon 
asymmetry due to the washout effects. A more reliable result is obtained if 
one uses the empirical fit [8] 

lo Sio d = iogiol 1 - + min{d!, d 2 , d 3 } , 

di = 0.8 log 10 k - 0.7 + 0.05 log 10 M x 10 , (31) 
d 2 = -1.2-0.05 \og w Ml° , 

d 3 = -(3.8 + log 10 Ml°)(\og w k - 1) - (5.4 - 0.67 \og w (Mi/GeV)) 2 - 1.5 , 

where £ is defined in Eq.(26), M\° = Mi/10 10 GeV and the function min 
evaluates to the smallest of the quantities dj. This fit reproduces considerably 
better the exact solution [6] of the Boltzmann equations for a wider range of 
K and Mi. We notice however that for Mi < 10 8 GeV and k ^> 1 the approx- 
imation (30) tends to give asymptotically a better result than the empirical 
fit in Eq. (31) [27]. 

We shall therefore use a combined fit in our analysis in order to estimate the 
washout effects. If Mi > 10 8 GeV or Mi < 10 8 GeV with « < 1, we calculate 

4 Another variable commonly used in the literature [6,22,8] is the mass parameter 
mi = (Mj,M D )ii/Mi ~ 1.1 x 1(T 3 k eV. 



10 



d using Eqs. (31). On the other hand, if Mi < 10 8 GeV and k > 1 we make use 
of the approximation (30). This will allow us to obtain a simple and reliable 
result for the solution of the Boltzmann equations without resorting to the 
full numerical solution of these equations. 



3.1 Analytic approach 



In this section we present a simple analytic approach that will allow us to ob- 
tain upper bounds for the baryon asymmetry in the present framework for the 
different solar neutrino solutions and assuming different patterns for the light 
neutrino spectrum. We shall divide our analysis in two parts. First we consider 
the case of large solar mixing, which includes four solar solutions: LMA, LOW, 
VO and Just-So 2 solar solutions. Secondly, we discuss the case of small solar 
mixing, i.e. the SMA solar solution. We shall consider for each case two possi- 
ble patterns for the light neutrino spectrum: hierarchical, i.e. mi <C rri2 <C m 3 
and inverted-hierarchical, mi ~ m 2 > m 3 , with the mass squared differences 
corresponding to the observed hierarchies Am 2 , 3> Am|. The analysis of an- 
other possible pattern, namely the case when the light neutrino masses are 
almost degenerate, mi ~ m 2 ~ m 3 , is more subtle and crucially depends on 
the inclusion of the CP- violating phases in the leptonic mixing matrix. There- 
fore we shall not discuss it here. However, it turns out from our full numerical 
study of Section 3.2 that in the latter case the produced baryon asymmetry 
is highly suppressed for all the solar solutions. 

The PMNS mixing matrix U v can be parametrized in the standard form [12] 



where c. 



C12C13 

— S12C23 

S12S23 - Ci2C23Si 3 e 



Cl2"S23 s 13 e 
15 



IS 



S12C13 



C12C23 ~~ Sl2S23 s 13 e 



iS 



-C12S23 ~~ s 12C 2 3Sl3e 



iS 



si 3 e 

•S23 c 13 
C23C13 



\ 



P 



(32) 



cos 6{j , Sij 



sinOij and P = diag (1, e ta , e t/3 ); 5 is a Dirac-type 
phase (analogous to that of the quark sector) and a, j3 are two physical phases 
associated with the Majorana character of neutrinos. 



We shall assume for the Dirac neutrino mass spectrum the SO(10)-motivated 
hierarchy given in Eq. (17), i.e. M D ~ M u . Moreover, since for the up-quark 



masses at GUT scale the hierarchy m v 
(e ~ 3 x 1CT 3 ), we shall write 



m. 



m t ~ e : e : 1 is verified 



d D = midiag(e 2 ,e, 1) . 



(33) 



Finally, in order to simplify our analytical discussion, we neglect the possible 
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misalignment between the charged-lepton and Dirac neutrino mass matrices. 
In other words, we assume Vl — 1 in Eqs. (16) and (19). As discussed in 
Section 2 such an approximation is reasonable in the context of SO(10), where 
one expects V L to be of the order of the CKM matrix. A more refined study 
which includes this effect will be given in the next section when we present 
our full numerical discussion. 

Eqs. (32) and (33) together with the light neutrino mass spectrum are therefore 
the only input parameters necessary in our further analysis. 



Case I: Large mixing (LMA, LOW, VO, Just-So 2 ) 



We consider maximal mixing in the 2-3 sector of the leptonic mixing matrix, 
i.e 023 = 7r/4, but keep the solar angle 9u as a free parameter in order to 
account for deviations from maximal mixing in the 1-2 sector. Since our goal 
is to obtain an upper bound for the baryon asymmetry, for the moment we 
can neglect any CP-violating phase, i.e. assume 5 — a — f3 — 0, and max- 
imize the asymmetry by simply replacing the imaginary parts in the lepton 
asymmetry (21) with their corresponding absolute values. Finally, since the 
mixing angle 9 13 is constrained by reactor neutrino experiments to be small, 
U e3 = | sin6> 13 | < 0.2 [28], the PMNS mixing matrix (32) can be approximately 
written in the form 



( 



c 



S Q 

75 ( c o - s Q U e3 ) 



U e ; 
1 

V2 



\ 



{ 75^0 - c & U e3 ) -^(c + s Q U e3 ) ^ J 



(34) 



where s & = sin 9u, c = cos Ou- The effective neutrino mass matrix given in 
terms of the light neutrino masses is easily obtained from Eq. (18). 

Now we introduce the symmetric matrix 

2 4 
711 ( 

M' = -^d^M^ , A = m lC 2 & + m 2 s 2 + m 3 U 2 e3 , (35) 
which can be expressed as 



M' 



1 


pe 


r e 2 


pe 


qe 2 


s e 3 


r e 2 


se 3 


te 4 



(36) 
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with the coefficients 
1 



V 

Q 

r 

s = 



V2A 
1 



[-mi (s + c e U e3 ) + m 2 s (c - s Q U e3 ) + m 3 U e3 ] , 



^[mi(«0 + c & U e3 ) 2 + m 2 (c - s G U e3 ) 2 + m 3 ] , 



1 



-[mic (s - c Q U e3 ) - m 2 s Q (c Q + s G U e3 ) + m 3 U e3 ] , 



V2A 1 

— [m^clU^ - s|) - m 2 (c| - s 2 Q U 2 3 ) + m 3 ] , 
* = JA^ mi ^ S& ~ c&Ue ^ 2 + m2 ( c © + ^o^) 2 + m 3 ] . 



(37) 



The eigenvalues of M' are easily obtained in leading order of e. We find 



Ai ~ 1 + pV , A 2 ~ (q - v ) e , A 3 ~ 



s (2pr — s) + q(t — r 2 ) — p 2 1 



q — p A 



(38) 



According to Eqs. (19) and (35), these eigenvalues are nothing but the inverse 
of the right-handed Majorana masses with the proper normalization coeffi- 
cient. Thus from Eqs. (35)-(38) we obtain 



M, = — 



A- 1 , 2 = 1,2,3. 



(39) 



Next, we can find the unitary matrix Ur which diagonalizes M' and determines 
the structure of the Dirac neutrino mass matrix through the relations (16) and 
(19), i.e. 



M D = d D U R . 



In leading order of e we obtain 



( l-^e 2 



U, 



R 



pe 



V 



pe 



r e 



1 _1_ P 2 t 2 , (pr-s) 2 9 
-1 + — e + 2(q-p^ 6 



pr—s 



ps-qr 2 

pr—s 

1 t 

5-P 

1 _ (pr-s) 2 9 
1 2( g - p 2)2 £ / 



(40) 



(41) 



It is now straightforward to estimate the baryon asymmetry. From Eqs. (23), 
(24), (33) and (39)-(41) we find 



Yn = 



3 dp 2 (q-p 2 -2r 2 )+2prs + r 2 (t-r 2 ) m 



t A 



327T g* 



1 + p 2 + r 2 



(42) 
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We can also write down the parameter k which together with the mass M\ 
controls the washout effects. Using Eq. (29) we find 

_ M P A(l +P 2 + r 2 ) 

1.7x8^2/0- • ^ 



As mentioned before we are particularly interested in two limiting cases, 
namely, when the light neutrino spectrum is hierarchical or inverted hier- 
archical. Below we shall consider each one of these situations. 



( i ) Hierarchical spectrum 

In this case m\ <C m 2 — \j 'Am| <C m 3 ~ ^ Am 2 . an d the coefficients in the 
matrix (36) are approximately given by 



m 2 s 2& + 2m 3 U e3 m 3 
V ' ~ ~ ~ 



2V2{m 2 s 2 Q + m 3 U 2 3 ) ' " 2(m 2 s| + m 3 U 2 3 ) ' 

-m 2 s 2Q + 2m 3 U e3 
2V2(m 2 s% + m 3 U 2 3 ) ' 



r ~ -^2g 20 + 2m 3 [/ e3 ^ (44) 



with s 20 = sin 26*i 2 . 

The right-handed Majorana neutrino masses read as 



Mi 



2 4 2 

mfe mf. 



2m t 2 e 2 , 2 rr2 x 2m 2 (4 V /A^| + V /A^^ 3 ) 
(m 2 4 + m 3 C/ e 2 3 ) - 



m 2 m 3 s 2 Q V > s% s jAm%Am 2 a 

M 3 ~ ^f® . (45) 
zmi 



We notice that the requirement M 3 < Mp implies a lower bound on the 
lightest neutrino mass: 

TTl 2 S 2 

m 1 > ~ 2 x lO" 7 eV , (46) 

for a typical value of m t ~ 100 GeV at GUT scale and a large solar mixing 
angle 9\ 2 ~ 7r/4. Obviously, deviations from maximal solar mixing will slightly 
modify this bound. 

We can also set an upper bound for the lightest right-handed Majorana mass 
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Mi. Indeed, from the first equation in (45) and taking U e3 = we find 



2 



Mi < - " . (47) 

With s and Amg given by the lower bounds in the case of the Just-So 2 solar 
(cf. Table 1) we find Mi < 1.3 x 10 9 GeV. We notice that this value is consistent 
with the requirement Mi < Trh, where Trh is the reheating temperature 
after inflation, which in turn is constrained to be below 10 8 — 10 10 GeV from 
considerations of the gravitino problem [22]. 

Moreover, for the mass ratios Mi/ Mi , % = 2, 3 we have 

M ^ mlr & sl Mi ^ 2m>i ^ Ml ( 48 ) 

M 2 2ml{r & s% + U 2 3 ) 2 ' M 3 m 2 s 2 Q Jkmi(r e s 2 Q + U 2 3 ) M 2 ' 

where 



/v x 1 2 

I Am 2 



r©=hr-f • (49) 



a , 



Let us also write down the approximate expression for the baryon asymmetry: 
3d m 2 ml 1 s 2 U 2 3 m 2 4 



i.e. 



32vr^ (m 2 s| + m 3 U 2 3 ) 2 {m 2 2 s% + m|t/ 2 3 ) v 



3 dm 2 u r & s 2 Q U 2 3 



(50) 



B ~ 32vr ^ V 2 (r 4 + f/ 2 3 ) 2 (r 2 4 + U 2 3 ) ' l ° J 

As a function of £/ e 3 , Yb reaches its maximum value when 

q 2 3/2 

U 2 e3 = ^f- • (52) 
Substituting this value into Eq. (51) we obtain the upper bound 



32ng* v 2 r s| 9 ' 2 X 10 r Q s 2 Q ' 



It is then clear that the larger the mass squared difference Am Q is, the higher is 
the suppression in the CP asymmetry e^. Thus, while the Just-So 2 vacuum 
oscillation solution is the most favoured in this framework, the LMA solar 
solution turns out to be highly disfavoured. 
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Table 1 

Constraints on neutrino masses and mixing angles coming from global analyses of 
the solar, atmospheric and reactor neutrino data [29,30]. The numbers in square 
brackets correspond to the best-fit values. 







Atmospheric 


and reactor neutrinos 








Am 2 (eV 2 ) 


tan 2 9 a 




u e3 




(1.4- 


- 6.1) [3.1] x 10~ 3 


(0.4 - 3.1) [1.4] 




< 0.2 


Solar neutrinos 






Am 2 (cV 2 ) 


tan 2 Q 


r e = 


(Am 2 /Am 2 ) 1/2 


LMA 


(2- 


30) [4.5] x 10~ 5 


(0.3 - 1) [0.4] 


(0.6- 


4.6) [1.2] x 10" 1 


SMA 


(4- 


10) [4.7] x 10~ 6 


(2-8) [4] x 10~ 4 


(2.6 - 


8.4) [3.9] x 10~ 2 


LOW 


(0.3 


- 2) [1.0] x 10~ 7 


(0.5 - 1.1) [0.7] 


(0.2 - 


1.2) [0.6] x 10~ 2 


VO 


(4- 


6) [4.6] x 10- 10 


(1.5-4)[2.4] 


(2.6- 


6.5) [3.8] x 10" 4 


Just-So 2 


(5- 


8) [5.5] x 10~ 12 


(0.5 - 2) [0.7] 


(2.9 - 


7.6) [4.2] x 10~ 5 



In Table 2 we give the upper bounds for the different solar solutions obtained 
using the present constraints on neutrino masses and mixing angles coming 
from global analyses of the solar, atmospheric and reactor neutrino data (cf. 
Table 1). We see that even if one neglects the washout effects (i.e. if one 
assumes d — 1), the LOW and LMA solutions have an upper bound which is 
below the lower bound of the observed baryon asymmetry. On the other hand, 
the VO and Just-So 2 solar solutions have an upper bound which lies inside 
the allowed range for the asymmetry. 

A more realistic bound can be obtained by including the dilution effects. From 
Eq. (43) we find 

M P m\s% + m\Ul z ^ 9.1 x 10 2 s 2 Q Am 2 & + Am 2 t/ e 2 3 
1.7 x 87rv 2 y /g:m 2 sl + m 3 U 2 3 ~ 1 eV s 2 G ^Am 2 G + y / Am 2 f/ 2 3 ' 

(54) 

Here we can distinguish two different regimes for the parameter k depending 
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Table 2 

Upper bounds for the clean (d = 1) and net (d / 1) baryon asymmetries in the case 
of hierarchical neutrinos. The bounds are obtained from Eq. (53), using the solar, 
atmospheric and reactor neutrino data of Table 1. 



Solution 


LMA 


LOW 


VO 


Just-So 2 


SMA 


Yb (clean) 
Y B (net) 


3.1 x 1(T 13 
1.4 x 10~ 14 


8.7 x 1(T 12 
5.1 x 1(T 13 


4.3 x irr 11 

2.5 x 1(T 12 


3.2 x irr 10 
2.0 x irr 11 


3.3 x 1(T 10 
1.7 x HT 11 



on the value of U e3 : 



K ~ 9 .lx 10 2 L^ for U 2 e3 « r|4 « r 4 , 



/Am 2 

k~9.1 x lO 2 ^ " for C/ e 2 3 > r 4 > 44 . (55) 



Thus, for large values of ?7 e 3 the parameter k is controlled by the mass squared 
difference of the atmospheric neutrinos and for all the large mixing solar solu- 
tions it lies in the range 35 < k < 70. However, to establish an upper bound 
on the baryon asymmetry we use the value of U e3 in Eq. (52) which maximizes 
Yb- Substituting it into Eq. (54) we have 



2 (Am 2 Am 2 ) 1/4 
6.4 x 10 2 - 2 ®± — ? 



k — 6.4 x 10 2 - 2 ^ — , (56) 

ev 

and according to Eq. (45), the mass of the lightest right-handed Majorana 
neutrino is approximately given by 

m 

Ml c -^^= . (57) 

Using these expressions to calculate the dilution factor d we obtain the upper 
bounds given in Table 2. These bounds are of course more stringent than the 
ones previously obtained neglecting the washout effects. 

In Fig. 1 we plot the asymmetry Y B as a function of U e3 for the mass of the 
lightest neutrino m x = 10~ 6 eV. Fig. la corresponds to the clean asymmetry, 
i.e. neglecting the dilution effects, while in Fig. lb the net baryon asymmetry 
(after including the washout effects) is plotted. The parameter k as defined 
by Eq. (29) is given in Fig. lc. We notice that for very small values of U e % 
this parameter is proportional to the mass squared difference of the solar 
neutrinos, k oc (Am^) 1 / 2 , and as U e3 increases, it tends to a common value 
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Fig. 1. The baryon asymmetry Yb as a function of U e 3 for the different solar solu- 
tions. Fig. (a) corresponds to the clean asymmetry, while in Fig. (b) the net baryon 
asymmetry, which includes the washout effects, is given. The parameter k defined 
in Eq. (29) and the mass of the lightest right-handed Majorana neutrino M\ are 
plotted in Figs, (c) and (d), respectively. The curves are given for the mass of the 
lightest neutrino mi = 10 -6 eV. The values for the neutrino oscillation parameters 
are taken from Table 1. 



for all the solar solutions, which is determined by the mass squared difference 
of the atmospheric neutrinos, i.e. k oc (Am 2 ) 1 / 2 (see Eqs. (55)). A similar 
behaviour is observed for the lightest Majorana neutrino mass M\ as can be 
seen from Fig. Id. In the latter case, M\ oc (Am!,) -1 / 2 for small values of U e3 
and Mi oc (Am 2 ) -1 / 2 for large values of U e3 (see also Eqs. (45)). 



At this point one may wonder whether the inclusion of a misalignment between 
the charged-lepton and Dirac neutrino mass matrices could change our con- 
clusions. Since in SO(10)-motivated scenarios such a misalignment is expected 
to be proportional to the CKM quark mixing matrix, let us then assume for 
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Fig. 2. The same plots as in Fig. 1 but for A = 0.15 and a Dirac phase 5 = tt. 
the unitary matrix V L in Eq. (9) the following CKM-type real matrix 



1 - A 2 /2 A 
-A 1 - A 2 /2 A\ 2 
-AX 2 1 



(58) 



with A ~ 1 and A of the order of the Cabibbo angle, i.e. A < 0.2. 

Following the same analytic approach we can obtain upper bounds for the 
baryon asymmetry. Since the relevant formulae are more cumbersome in this 
case, it is more illustrative to present the results in a simple plot, similar to 
that of Fig. 1. The results are presented in Fig. 2. From this figure it is seen 
that the upper bounds previously found for Yg are essentially unaltered and 
that our previous conclusions remain valid in this case. To maximize the effect 
we have taken the Dirac phase 5 in the PMNS neutrino mixing matrix to be 7r, 
as suggested from our full numerical study (see next section). The resonance 
behaviour in the curves at U e s ~ 0.1 is associated with our particular choice 
for the values of S and A. 



( ii ) Inverted-hierarchical spectrum 
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In this case m 3 «mi~m 2 ~ ^Am^ and the coefficients in the matrix (36) 
are approximately given by the expressions: 

_ (m 2 - m 1 )s 2Q _ U e3 
P ~ 2 v / 2(m 1 c| +m 2 s|) ' 
^ mis| + m 2 c _ (m 2 - mi)s 2Q Z7 e3 
2(mic|+m 2 s|) 4(mic| + m 2 s|) 
(m 2 - mi)s 2 ^3 



2 v / 2(m 1 c| + m 2 s|) ' 
_ mis| + m 2 c| 1 2 

2(m lC | + m 2 4) ^2 e3 ' 
_ mis| + m 2 c| (m 2 - mi)s 20 Z7 e3 



with 



2(miCg + m 2 s|) 4(mic| + m 2 s|) 
m 2 — mi 1 Airig r 



mic| + m 2 s| 2 Am^ 2 
The right-handed Majorana neutrino masses are given in this limit by 



(60) 



2 4 2 

rare mf, 



mi 4 + m 24 jAm 



Mi 



M 2 ~ 2m 2 e 2 f ^ + ^ ~ -%lL 
V m 2 mi/ ./Am! 



M 3 ~ -^j mi+m2 ~ . (61) 
4m 3 mi4 + m 2 4 2m 3 

Once again requiring M 3 < M P and taking m t ~ 100 GeV at GUT scale, we 
find the following lower bound on the lightest neutrino m 3 , 

2 

777 

ms ~ 2M~ " 4 X 10 " 7 eV • (62) 



The generated baryon asymmetry can be approximated in this case by the 
following expression 

Y 3 d m x m 2 {m 1 - m 2 ) 2 s 2 2& + 2m 3 (mf + ml])Z7 2 3 m\ ^ 
B 128ir g* (m\c% + m 2 s%)(mic% + m 2 s%) 2 v 2 



or equivalently, 



3 d m 2 . 

b — 

5127?^ v 2 



2 4 , 16m 3 f/ 2 3 

A 20'O "I" 



(64) 
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Taking for instance m 3 ~ yAm| <C — y Ami an d the maximum value 
of U e 3 ~ 0.2 it is easy to see that the dominant contribution to the asymmetry 
comes from the second term in the square brackets. Neglecting the washout 
effects (d = 1) we find the following upper bound: 

Yb < ^-^eVo - 3.7 x 10- 16 r o . (65) 



Thus we conclude that the baryon asymmetry is highly suppressed for all the 
large mixing solar solutions and for an inverted hierarchical neutrino spectrum, 
even without taking into account the washout effects. It is however worthwhile 
to have an idea of these effects. The parameter k reads as 



Mr 



m\ + m\ 



1.7 x IQirv 2 , 



M, 



fglmic% + m 2 s© 



1.7 x Snv 2 , 



(66) 



Therefore, 



k — 9.1 x 10' 



v' 



Am 2 a 



eV 



~ 50 . 



(67) 



Since from Eq. (61) we have Mi ~ 1.8 x 10 4 GeV, then from Eqs.(31) we find 
10" 2 . 



Case II: Small mixing (SMA) 



Let us now consider the small mixing solar solution. Most of the relevant 
formulas can be easily obtained from the previously discussed large mixing 
case by just letting c — > 1 and assuming s C 1. We must however proceed 
with care since in this case there are two small parameters competing against 
each other, namely, the small solar mixing angle 9 12 and [/ e3 . In this case the 
leptonic mixing matrix U v can be approximated as 



{ 73(^0 - u e3 ) 



72(1-^3) 75 
-4^(1 + ^3) js) 



(68) 



Next we follow the same steps performed in the large mixing case. We first 



21 



define the new matrix 

2 4 

M' = T ^-d D 1 M u d D 1 , A = m l + m 2 s% + m 3 U 2 e3 , (69) 

which can be written in the form (36) with the coefficients p, q, r, s, t as given in 
(37). The eigenvalues and eigenvectors of M' are given by the same expressions 
in Eqs. (38) and (41), respectively. The right-handed Majorana masses read 
as in Eq. (39). Finally, the baryon asymmetry will be given by Eq. (42) with 
the parameter k defined in Eq. (43). 

Let us now consider the limiting cases of hierarchical and inverted-hierarchical 
neutrino mass spectrum. 



( i ) Hierarchical spectrum 

In the limit mi <C m 2 <C m 3 we find 

m 2 s Q + m 3 U e3 m 3 



V 



y/2 (mi + m 2 s| + m 3 U 2 3 ) ' ' ~~ 2(m l + m 2 s% + m 3 U 2 3 ) 

-m 2 s Q + m 3 U e3 ^ 

y/2 (mi + m 2 s| + m 3 U 2 3 ) 



Moreover, 
Mi 



m?e 4 m 2 



mi + m 2 s% + m 3 U 2 e3 mi + ^Amf s| + ^Am 2 a U ( 



^ 2m 2 t e 2 mi + m 2 s% + m 3 U 2 3 ^ 2m 2 c m 1 + ^ Am% s% + ^/Am 2 a U 2 3 
m 3 mi + m 2 (s Q - U e3 ) 2 ~ f&~2 mi + J Am% (s - U e3 ) 2 



M 3 ^(± + (SQ ~ U «A a ^ -J— + (S0 = ^ 1 (71) 
2 \^m 2 mi y 2 ^Am| m i 



From the requirement that M 3 be smaller than the Planck mass we find the 
following lower bound for the mass of the lightest neutrino: 



m? JAm%(s Q - U e3 ) 2 

mi > <V ^- ~ 4.1 x 10- 7 (s - f/ e3 ) 2 eV . (72) 

2MpJAm% - m? 



Using for instance the minimum value s = 0.014 and taking U e3 = one has 
m!>8x 10~ n eV. The above bound is of course sensitive to the value of U e3 
and it gets weaker as U e3 approaches s . 
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Let us now estimate the baryon asymmetry. From Eqs. (42) and (70) we find 

mg ( m i + m 2s|) U 2 3 e 4 



3 d m? 

327T g* v 2 (mi + m 2 s| + m 3 U 2 3 ) 2 (mf + m\s% + m|f/ e 2 3 ) 

3 dml (Amir/ 2 (m 1 + sl^A^)U 2 e3 _ 

-2 

e3j 



32tt ^ w 2 ( mi + s 2 Q J Am 2 & + \MmM) 2 (m? + s|Am| + Am 2 ?/, 



(73) 



We can also compute the parameter k relevant for the dilution effects. From 
the definition (43) we find 



K 



m\ + m\s 2 & + m\U { 



e3 



1.7 x % r KV 2 J~glmx + m 2 S0 + m 3 U 2 3 



(74) 



As a function of ?7 e 3 the expression (73) reaches its maximum when 

(m! + Ami) 1 / 2 (m 2 + s^Aml) 1 / 2 



17. 



e3 



v/2(Am 2 ) 3 / 4 



For m 1 3> s \/ Ami ^ s o v Am® we have 



r/2 _ 

63 " V2 



3/2 



mi 



Am 2 



and therefore 



3 dml\Am 



32ir g* v 2 mi 
For the parameter k we have in turn 

-\2 



~ 7.2 x l(T ib d 



eV 



mi 



6.4x10 1/9/. 9 /777i 

« ~ i ^ m} /2 (Am 2 ) 1/4 ~ 1.8 x lO 2 ^/-^ . 
Finally, for the lightest right-handed Majorana mass we obtain 

,2 reV - 



777 

Mi ~ — * ~ 10 3 GeV 
mi 



mi 



(75) 



(76) 



(77) 



(78) 



(79) 



It remains to estimate the dilution factor d. Assuming e.g. mi = 10 eV, 



so that the relation m x ^> s Q \J Ami lfi satisfied, Eqs. (78) and (79) imply 
k ~ 1.8, Mi = 10 7 GeV, and therefore d ~ 0.14. Thus, from Eq. (77) we 
obtain the upper bound 



Yb < 10 



-12 



(80) 
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which is one order of magnitude below the required value for the observed 
asymmetry. 



On the other hand, if m 1 <C s|^/Am| <C s Q Am^ then the maximum value 
in Eq. (73) is obtained when 

Ul = ^£ • (81) 

In this case, the expression for the upper bound of the baryon asymmetry is 
the same as the one previously found for the large mixing solar solutions (see 
Eq. (53)): 

^<7j--^^T^1.8xl(n 9 d. (82) 



The dilution factor turns out to be crucial in obtaining a reliable estimate for 
Yg. In this limiting case, the parameter k is also given by the same expression 
of the large mixing solutions, i.e. by Eq. (56). We find k ~ 8. For the lightest 
Major ana mass we have in turn, 

2 

M 1 ~ ~ 2.5 x 10 9 GeV , (83) 

4 v Am o 

and the mass ratios M\/M 2 and Mi/M 3 are given by the approximate expres- 
sions: 

M 2 ~ 2 \mj r s| ' M 3 ~ \m t ) 4 ' 



The above values for /t and Mi imply then a dilution factor d ~ 8 x 10 3 . 
Substituting this value into Eq. (82) we conclude that for the SMA solution 

Y B < 1.4 x 10- 11 , (85) 

a bound which is within the allowed range for the baryon asymmetry. 



( ii ) Inverted-hierarchical spectrum 

Let us now consider an inverted-hierarchical neutrino spectrum, i.e. mi ~ 
m 2 ~ J Ami ^ m 3- I 11 thi s case the coefficients of the matrix (36) can be 
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obtained from Eqs. (59) by setting c = 1 and assuming s <C 1. We find 

(m 2 -mi)% Ue3 s Q rl U e3 m 2 1 

p ~ = -j= ~ — ^ -j= , g ~ — s ~ t ~ 



V^mi V2 2^2 v 7 ^ ' ' 2mi 2 ' 

(86) 



(m 2 - mi)% C/ e3 s r 2 [/ e3 



v^mi ^2 2^ 72 ' 



The right-handed Majorana neutrino masses will be given by the same ex- 
pressions in Eqs. (61) and, in particular, the lower bound (62) for the lightest 
neutrino should be verified as well. Finally, the approximated expression 
for the baryon asymmetry is obtained from Eq. (64): 

Yb - ;J--^^t=| < 3.7 x l(T 16 r < 3.1 x 1CT 17 , (87) 
327T 5 * v 2 J Am* 



for d = l,m 3 ~ yArriQ and U e % = 0.2. Thus we conclude that for the SMA 
solution and an inverted-hierarchical neutrino spectrum the generated baryon 
asymmetry is highly suppressed. 



3.2 Numerical Analysis 



From the simple analysis performed in the last section we have concluded that 
in the SM framework with the simplest SO(10)-motivated hierarchy (17) for 
the Dirac neutrino mass spectrum, only the SMA and Just-So 2 solutions of 
the solar neutrino problem are compatible with the required value for the 
baryon asymmetry of the Universe. In this section we will perform a full 
numerical computation of the baryon asymmetry including not only all the 
CP- violating phases which appear in the PMNS leptonic mixing matrix U v 
defined in Eq. (32), but also the possible misalignment which may result from 
the process of diagonalization of the charged-lepton and Dirac neutrino mass 
matrices. This misalignment is characterized by the matrix Vj, in Eq. (16), for 
which we shall assume the CKM structure of Eq. (58). 

We then proceed as follows. The PMNS mixing angles as well as the mass 
squared differences Am Q and Am 2 , are allowed to vary in the ranges indi- 
cated in Table 1. First we randomly fix a set of values for the PMNS mixing 
angles and the CP-violating phases a, (3 and 5. The latter are randomly 
chosen in the interval from to 2ir. In order to compute the effective neu- 
trino mass matrix we take as an input the mass of the lightest neutrino, mi, 
which varies from 10~ 10 to 1 eV for the SMA solution and from 10~ 7 to 1 
eV for all the large mixing solutions^ This means that unlike the analytic 

5 The lower bounds for the mass of the lightest neutrino are taken in agreement 
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(a) LMA 
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m, ( eV ) m, ( eV ) 



Fig. 3. The baryon asymmetry as a function of the mass mi of the lightest neutrino 
for the LMA, SMA, LOW and VO solar solutions. The black-dotted area corre- 
sponds to sets of input parameters which yield an asymmetry within the allowed 
range given in Eq. (2) and delimited by the dashed lines. 

study, the full numerical analysis includes also the case of almost degener- 
ate neutrinos. This is possible because we are taking into consideration all 
the phases that may be present and, consequently, possible cancellations. As- 
suming mi < m 2 < m 3 , the values of wi2 and m 3 are determined from the 
expressions m 2 = \Jm\ + Am| , m 3 = \jm\ + Am| + Ami ■ The inverse hi- 
erarchical case is not discussed here since it leads to a very suppressed baryon 
asymmetry as already shown in the analytic discussion. The effective neutrino 
mass matrix M v in Eq. (18) is computed, and also the matrix M in Eq. (19) 
with the up-quark masses given at the GUT scale: m u ~ 1 MeV, m c ~ 0.3 GeV 
and m t ~ 100 GeV [31] and the matrix V L in Eq.(58) with A randomly cho- 
sen in the interval < A < 0.22. The matrix MM^ is then diagonalized to 
obtain Ur and the right-handed neutrino masses Mj. In the basis where Mg 
and Mr are diagonal the Dirac neutrino mass matrix Mo is determined by 
Eq. (16). Finally we compute the baryon asymmetry from Eqs.(21), (22), (25) 
and (27), where the dilution factor d is determined using a combined fit from 
the approximations given in Eqs. (30) and (31). 

with our previous analytic estimates, so that the requirement M3 < Mp is always 
satisfied. 
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71/2 % 331/2 2ti 0.05 0.1 0.15 0.2 



8 X 

Fig. 4. In Fig. (a), the baryon asymmetry as a function of mi in the case of the 
Just-So 2 solar solution. The black-dotted area corresponds to sets of input param- 
eters which give an asymmetry within the allowed range (dashed lines) as given 
in Eq. (2). In Fig. (b), the right-handed Majorana neutrino masses are plotted as 
functions of m\. From Fig. (c) we conclude that for relatively large values of U e s 
the Dirac phase 5 should be close to ir to produce a baryon asymmetry in the al- 
lowed range. The allowed region in the plane (A, U e3 ) is shown in Fig. (d). The 
gray-dotted area is obtained by slightly relaxing the lower bound on the baryon 
asymmetry, Yb > 1CT 11 . 

In Figs. 3 we present the results for the baryon asymmetry Yb as a function of 
mi for the LMA, SMA, LOW and VO solutions of the solar neutrino problem. 
The first immediate conclusion which can be drawn from these plots is that 
among these four solutions only the SMA solution is compatible with the 
experimental range for Yb with a mass for the lightest neutrino in the range 
from 2 x 10~ 10 to 4 x 1CT 6 eV. Moreover, the values of U e % should be small in 
this case, \U e3 \ < 1CT 2 . By comparing the results plotted in Figs. 3 with the 
upper bounds for Y B obtained in the last section (see Table 2), we find the 
analytic bounds in very good agreement with the exact numerical ones. 

Finally, in Figs. 4 we present the results for the only large mixing solar solu- 
tion that produces the required Yb, namely, the Just-So 2 vacuum oscillation 
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Fig. 5. The CP-violating Majorana phases a and (3 in the cases of no Dirac-type 
CP violation (Figs, (a) and (b)) and maximal Dirac-type CP violation (Fig. (c)) for 
the Just-So 2 vacuum oscillation solar solution. The black-dotted area in the plots 
refers to the values which lead to an acceptable baryon asymmetry Yb- 



solar solution. The asymmetry as a function of the lightest neutrino mass mi 
is plotted in Fig. 4a. The black-dotted area corresponds to sets of input pa- 
rameters which give an asymmetry within the allowed range delimited by the 
dashed lines (see also Eq. (2)). In Fig. 4b, the right-handed Majorana neutrino 
masses are plotted as functions of m\. As before, the lower bound on m\ is 
determined by requiring the mass of the heaviest Majorana neutrino M3 to be 
below the Planck scale Mp. From Fig. 4c we conclude that for |C/ e3 | > 10~ 2 
the Dirac phase S should be close to n to produce a baryon asymmetry in the 
allowed range. This in turn corresponds to larger values of the misalignment 
parameter A in the matrix (58) as can be seen from Fig. 4d. 

From Figs. 3 and 4 it is also clear that an acceptable baryon asymmetry 
requires very small values of the mass of the lightest neutrino, mi C 1 eV. 
This in particular implies that, in the present framework, an almost degenerate 
neutrino spectrum is excluded. 

To establish a relationship between CP violation at high and low energies is 
not an easy task in the sense that a general analytic treatment with the full set 
of parameters would be difficult to perform. Nevertheless, one can investigate 
whether there is any correlation between the CP-violating phases a, (3 and S in 
a plausible leptogenesis scenario. In Fig. 5 we present the results of a random 
analysis in the cases of no Dirac-type CP violation (5 = 0, ir) and maximal 
Dirac-type CP violation (5 = n/2) for the Just-So 2 vacuum oscillation solar 
solution, which is the only large mixing solar solution compatible with lep- 
togenesis in the simplest SO(10)-motivated framework considered here. The 
black-dotted area in the plots refers to the values of a and j3 which lead to an 
acceptable baryon asymmetry Y#. It is seen from Figs. 5a and 5b that one can 
obtain the right magnitude for Y B with a single nonvanishing Majorana phase. 
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This means that the possible absence of CP violation effects in neutrino os- 
cillations in future experiments may not a priori discard leptogenesis as the 
mechanism responsible for the generation of the BAU. In the case of maximal 
CP violation in Fig. 5c, the phases a and (3 are slightly less constrained than 
in the previous case. 

We have seen that in the present framework and assuming a mass hierarchy for 
the Dirac neutrinos like the one of up-quarks (cf. Eq. (17)), it is not possible 
to reconcile the LMA, LOW and VO solar solutions with the leptogenesis 
scenario. We may then ask ourselves what type of hierarchies should the Dirac 
neutrino masses satisfy for these solutions to be realized and leptogenesis to 
be viable. To answer this question let us now relax our previous assumption 
(17) and write the most general form for the Dirac neutrino mass spectrum: 

d D = m D3 diag(ei, e 2 , 1) , (88) 

with £j = mr)i/mD3, i — 1,2 and < e± < £2 < 1- 

For a given scale of the heaviest Dirac neutrino moz , we can vary the param- 
eters (£1,62) i n the above interval and look for the allowed regions where the 
cosmological baryon asymmetry is within the experimental range. To illustrate 
the results, we shall consider only the presently most favoured solar solution, 
i.e. the LMA solar solution. The LOW and VO solutions can of course be 
analyzed in a similar manner. The results are plotted in Fig. 6 for two typical 
scales, namely, mm of the order of the bottom-quark mass m& ~ 1 GeV and 
around the top-quark mass m t ~ 100 GeV at GUT scale. It is then clear from 
the figure that one may obtain the correct baryon asymmetry, provided one 
uses a hierarchy for the eigenvalues of corresponding to £j lying inside 
the allowed ranges indicated in Fig. 6. However, note that for Md oc M u or 
Mo oc Md one is not able to have the LMA solar neutrino solution and yield 
the correct baryon asymmetry. 

We emphasize that this analysis is quite general and does not rely on any spe- 
cific texture for Dirac matrix Mp or the Majorana matrix Mr, except for the 
fact that in the basis where the charged leptons and right-handed Majorana 
neutrinos are diagonal, we assume a CKM-type misalignment between M^ and 
M D . 



4 Implications for low energy CP violation and neutrinoless double 
beta decay 

It is clear from the analysis presented in the previous sections that CP viola- 
tion is a crucial ingredient for the generation of the cosmological baryon asym- 
metry through leptogenesis. In our particular framework, the factor M D M D 
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Fig. 6. Allowed region for the Dirac neutrino mass ratios e% = mpi/mps in the case 
of the LMA solar solution. We assume m\ = 1CT 3 eV and two different scales for 
the heaviest Dirac neutrino, tuds ~ rrib = 1 GeV and mm ~ mt = 100 GeV. 

is of the form 

Ml M D = (U R d D V L ) (y£ d D U R ) = U R d 2 D U R . (89) 

As a result, it does not depend directly on Vl, yet the presence of Vl is felt 
in the determination of Ur through Eqs. (19) and also of the heavy Majorana 
masses (i=l,2,3), which in turn appear in the lepton asymmetry e^. In 
Eq. (89) with Ur written as U' P\ (in the notation of Eq. (14)) it is obvious 
that the three phases in Pg do not play any role for leptogenesis. Moreover, in 
the limit Vl — 1, we have 

M D = d D U R = P ( d D U' P 1 , (90) 

and in this case P^ can be rotated away from M D so that only three physical 
CP- violating phases remain in M D . In this case it is also possible to choose a 
WB where both M D and Mi are diagonal and real, hence 

M R = U* R d R Ui, (91) 

implying that all CP- violating effects could be generated at high energies with 
CP softly broken. 

One of the striking features of our numerical results is the fact that in the 
present framework leptogenesis favours a Dirac phase 5 in the PMNS mixing 
matrix U u very close to n for relatively large values of \U e3 \ > 10~ 2 . In the 
exact limit 5 = ir all CP- violating effects at low energies are due to the 
presence of the Majorana phases a and/or (5. Nevertheless, the determination 
of the nature of the phases relevant to the lepton asymmetry is not trivial 
and would require a careful analysis of the matrix R in Eq. (10). This is due 
to the fact that Ur and U v are related in a complicated manner even in the 
limit V L — 1 as can be seen from Eqs. (19). Another remarkable feature is 
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the fact that the Majorana phases a and (3 can only give rise to the necessary 
baryon asymmetry in those regions of the parameter space represented by the 
bands in Figs. 5, where one can see that although all possible values of a and 
P are allowed, they are not independent of each other. The type of pattern 
generated in the (a, /3)-plane remains very similar for different values of the 
Dirac phase 5. It is also noticeable that leptogenesis is viable with a single 
nonzero CP-violating phase of Majorana character even for 5 = 0, tt, whilst 
the Dirac phase alone is not sufficient. 

The three phases 5, a and (3 contained in the mixing matrix U v (for 3 gener- 
ations of light neutrinos) are precisely the ones related to CP- violating phe- 
nomena at low energies. It is well known that neutrino oscillations are only 
sensitive to Dirac-type phases since the combination of matrix elements of 
U v appearing in the oscillation probability is such that the Majorana phases 
cancel out. CP violation in oscillations will manifest itself in the difference of 
CP-conjugated oscillation probabilities, e.g. 

P{Ve -> ^) - Pipe -> P„) OC J , (92) 

with J a function of the PMNS matrix which can be expressed in terms of 
experimentally relevant quantities as 

J = isin2# sin2# a (l - |f/ e3 | 2 ) |C/ e3 | sin5. (93) 



In the quark sector of the SM the parameter J defined in terms of the Vqkm 
elements is known to be 0(A 6 ) ~ 10~ 5 . One might expect that for large solar 
mixing J would be much larger, yet the constraints imposed by leptogenesis 
in the context of SO(10) with the Dirac neutrino masses fixed by the up quark 
spectrum imply an upper bound on J of the order of 10~ 4 , which is outside 
the experimental reach of the next generation of neutrino experiments. We 
notice however that within the LMA solution, which is at present the most 
favoured after the recent SNO results, and allowing for the most general form 
for the Dirac neutrino spectrum as in Eq. (88), one can reach values of J 
of the order of 1CT 2 , thus rendering CP-violating effects visible in the near 
future. 

Neutrino oscillation experiments give no evidence about whether neutrinos are 
Dirac or Majorana particles. Processes which violate lepton number such as 
neutrinoless double beta decay ((/3/3)oi,-decay) of even-even nuclei, (A, Z) — > 
(A, Z + 2) + e~ + e~, would imply that neutrinos have Majorana character. 
The probability amplitudes of these processes are proportional to the so-called 
"effective Majorana mass parameter" 



\<m>\ 



m 1 U e 2 1 +m 2 U e l + m 3 U e 2 3 , (94) 
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(with the possibility of complex C/ ei ), or equivalently, 



\<m>\ 



mi 4 (1 - |£/ e3 | 2 ) + m 2 4 (1 - |f/ e3 | 2 ) e 2M + m 3 |f/ e3 | 2 e~ 2 * 5 e 2 ^ 



(95) 



Presently, the most stringent constraint on the parameter | < m > | comes from 
the 76 Ge Heidelberg- Moscow experiment [32], indicating that |<m>| <0.35 
eV (90 % C.L.). This value is still too high to provide some information about 
the neutrino spectrum. In fact, it can only rule out some models with quasi- 
degenerate Majorana neutrinos. Nevertheless, higher sensitivity is planned to 
be reached in new generation experiments. The NEM03 [33] and CUORE [34] 
experiments intend to achieve a sensitivity up to |<m>| ~ 0.1 eV while for 
the EXO [35] and GENIUS [36] experiments this value is about one order of 
magnitude below. 



A hierarchical spectrum with mi < m 2 — yAm 2 <C m 3 ~ y Am 2 , together 
with the CHOOZ bound on \U e3 \ would imply |<m>| < 10 -3 eV, far below 
the present bound. On the other hand, for an almost degenerate spectrum 
with m = mi ~ m 2 — m 3 , for m > 0.35 eV some cancellation between the 
terms in (95) is already required, while in the case of small solar mixing, the 
contribution comes mainly from the term proportional to m x . Finally, for an 
inverted-hierarchical spectrum where m 3 < mi ~ m 2 — \J Am 2 a , the maximal 

possible value is |<m>| ~ \J Am 2 a ~7x 10~ 2 eV [37]. 

We conclude this section with the following remark. It has been recently 
claimed [38] based on a refined analysis of the data of the Heidelberg- Moscow 
experiment that there is evidence up to 3.1a- for the observation of neutrino- 
less double beta decay with |<m>| = (0.05 - 0.84) eV (at 95% CL) with 
a best-fit value of 0.39 eV. If confirmed, this result would imply an almost 
degenerate or inverted hierarchy for the light neutrino masses in the case of 
three generations [39]. If this is the case, baryogenesis via leptogenesis based 
on the simplest SO(10) GUT scenario would not be a viable mechanism, since 
it is not possible to reconcile any of the solar solutions with the required 
cosmological baryon asymmetry, as it follows from our analysis. 



5 Conclusions 



We have studied leptogenesis in the framework of a simple extension of the 
SM where the right-handed neutrinos are added to the standard spectrum, 
thus leading to the seesaw mechanism. In order to restrict the number of free 
parameters, we have made use of some GUT-inspired relations for the quark 
and lepton mass matrices. We have shown that the latter relations, together 
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with the constraints from the low energy neutrino data, imply important re- 
strictions on the size of leptogenesis. In particular, we have pointed out that 
for the Just-So 2 and SMA solar solutions, one can generate sufficient BAU 
through leptogenesis even for the simplest SO(10) GUT. On the contrary, for 
the LMA, LOW and VO solar solutions, a different hierarchy for the Dirac 
neutrino masses is required in order to obtain a viable leptogenesis. 

We expect our analysis and conclusions to remain valid also in the super- 
symmetric version of the present framework. Although in this case new decay 
channels will enhance the generated CP asymmetry, these additional contri- 
butions tend to be compensated by the washout processes which are in general 
stronger than in the SM case [40]. 

A related and very important subject is the search for CP violation in the 
leptonic sector, at low energies. This is at present one of the great challenges 
in particle physics and it has recently received a great deal of attention. Ex- 
periments with superbeams [41] and neutrino beams from muon storage rings 
(neutrino factories) [42] have the potential to measure directly the Dirac phase 
S through CP and T asymmetries [43] or indirectly through oscillation prob- 
abilities which are themselves CP conserving but also depend on 5 . An alter- 
native method proposed and discussed recently [44] is to measure the area of 
unitarity triangles defined for the leptonic sector [45]. 

In this paper, we have investigated the possible connection between CP viola- 
tion at low energies and leptogenesis. In our SO(10) inspired framework, once 
masses and mixing for three light neutrinos are fixed, the baryon to entropy 
ratio Yb can be computed. It is well known that CP violation at low ener- 
gies depends on three physical phases appearing in the PMNS matrix. One of 
them, 5 , is a Dirac-type phase to which neutrino oscillations are sensitive, the 
other two are Majorana-type phases which are relevant to processes such as 
neutrinoless double beta decay. We have shown that the required cosmological 
baryon asymmetry can be produced with only one of the Majorana-type phases 
different from zero, while the same is not possible with only a nonvanishing 
Dirac-type phase. Concerning the important question of whether the strength 
of CP violation at low energies will be sufficient to be measured through neu- 
trino oscillations, we have verified that there are two possible scenarios: if 
one assumes minimal SO(10) then, the constraints imposed by leptogenesis 
render these effects too small to be seen in the next generation of neutrino 
experiments; on the other hand, if one allows for a more general spectrum 
of the Dirac neutrino mass matrix, then the strength of CP violation can be 
sufficient to be visible at low energy neutrino oscillation experiments. 

The leptogenesis scenario crucially depends not only on the mechanism that 
was responsible for populating the early Universe with right-handed neutrinos 
but also on the precise details of the reheating process after inflation. In this 
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paper we have only considered the conventional scenario, namely, the decays of 
right-handed Majorana neutrinos which are produced in thermal equilibrium 
processes. Other production mechanisms may be as well viable and even com- 
petitive [22]. Their realization will of course depend on the particular details 
of the inflationary epoch and the thermal evolution of our Universe. 
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